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Abs  t  rac  t 

The  ESFRIT  algorithm  has  been  shovn  to  be  an 
effective  solution  to  the  angle  of  arrival  estima¬ 
tion  problem  [l].  One  possibility  for  the  imple¬ 
mentation  of  ESPRIT  is  the  use  of  a  linear  array 
to  provide  for  the  sensor  pairs.  This  paper  dis¬ 
cusses  a  technique  for  compensation  of  mutual 
coupling  effects  between  array  elements.  Computer 
simulations  demonstrate  a  significant  improvement 
in  performance. 

Introduction 

Direction  finding,  which  involves  estimation 
of  the  angles  of  arrival  of  sources,  is  very  im¬ 
portant  in  many  sensor  systems  such  as  radar, 
sonar,  seismology,  etc.  Several  au'hors  have  ap¬ 
proached  the  problem  using  subspace  methods  |1,2]. 
However,  these  methods  have  not  taken  into  account 
effects  of  mutual  coupling  between  array  elements 
which  can  significantly  deteriorate  the 
ergensystems  underlying  the  solution  procedures. 
In  this  paper  we  deal  with  compensation  of  the 
mutual  coupling  effects  when  a  linear  array  con¬ 
sisting  of  m  sensors  is  used  in  conjunction  with 
the  ESPRIT  algorithm  (  1  ] .  The  method  of  moments 
[3,4]  is  used  to  obtain  the  matrix  of  mutuals  for 
each  sensor  pair.  A  transformation  matrix  is  de¬ 
veloped  which  processes  the  observed  data  so  as  to 
estimate  the  signals  that  would  have  resulted  had 
there  been  no  mutuals.  We  show  that  ideally  the 
effects  of  mutual  coupling  can  be  completely 
eliminated.  Computer  simulations  demonstrating  the 
improved  performance  are  presented. 

Mutual  Coupling 

Consider  a'  linear  array  of  m  dipoles 
uniformly  spaced  at  a  distance  D.  Each  dipole  is 
of  length  t  and  has  a  radius  r  satisfying  the  con¬ 
dition  r<<t.  A  load  is  attached  to  the  center  gap 
of  each  dipole.  Assume  there  are  d  narrowband  sig¬ 
nals  impinging  on  the  array  as  planar  wavefronts. 
The  voltages  induced  by  the  assumed  signals  on  the 
ioads  are  the  outputs  of  the  dipoles.  Induced  cur¬ 
ler,  ts  v!  11  appear  on  the  dipoles.  These  currents 
reradiate  and  generate  wenffered  fields.  The  scat¬ 
tered  fields  then  induce  currents  on  ,r.e  neighbor¬ 
ing  dipoles.  The  process  of  induction  and  reradia¬ 
tion  causes  mutual  coupling  between  the  dipoles. 
Using  single  sinusoidal  expansion  and  weighting 
functions  per  dipole,  the  method  of  moments  [3.4| 


is  employed  to  obtain  the  matrix  of  mutuals. 
Denote  the  current  distribution  in  the  array  of 
dipoles  by  J(z)  (assuming  longitudinal  distribu¬ 
tion  and  neglecting  all  other  distributions)  ar.d 
the  j-th  expansion  function  by  fj(z).  Then 
m 

J(Z)*E  I(j)f;(z)  (1) 

3=1 

where  I(j)  are  unknown  amplitudes  to  be 
determined.  At  a  point  fy,z)  in  the  Y-Z  plane,  the 
scattered  field  is  given  by 

m 

E(s>(y,z)=r  I(j)E<f>(y.z)  (2) 

j-1 

where  Ej(s)(y,z)  is  the  scattered  field  from  the 
jth  dipole. The  total  field  is 

E(y,z).E(inc>(y,z)  *  E(s)(y,z)  (3) 

where  E^lnc^  is  the  incident  field.  Let  Ez  be  the 
z-component  of  the  total  field.  A  generalized 
voltage  V(i)  induced  on  the  subsection  spanned  by 
the  function  fj(z)  can  be  defined  with  respect  to 
a  weighting  function  Vj(z)  as 

V(i)=F(Ez(y,z),Wj(z))  (4) 

where  F  is  bilinear  with  respect  to  Ez(y,z)  and 
w;(y,z).  Similarly,  we  define 

v(inc)(l)»F(Epnc>(y,z),wi(z)),  (5) 

V<s>(i)=F(E|s)(y,z),Vj(z)).  (6) 

Thus, 

V( ( )=V< inc>(i ).V<S)( i ), 

which,  for  metallic  scatterers.  becomes 

V(i)=v(i,lc>(i)*v(s)(i)=0  , 

V(inc)(i)=-V<s>(i)  .  (7) 

However , 


m 

V<S>U)  »F(  r  I(j)E<f>(v,z).Vj(z)) 

j’l 


1 


m 

=  r  i(j)F(E(?)(y,z)>wi(z)) 

J-i 


z'j  — F(E(|>(y,z),wi(z)).  (8) 

Then 

m 

V<s>(i)  =  E  -zlJ  I(j)  ;  i  =  1 , 2 . m.  (9) 

j-1 


In  matrix  notation  V^s)=-Z  I  where 


Y=H“'x*M.  (13) 

Let  G=H'^.  G  can  be  written  as 


'  811 

812 

•  •  •  8 1 ( m+ l ) 

821 

822 

•  •  ■  82 ( m+ 1 ) 

■  &(rti.l)l 

8( m+ 1 ) 2  • 

•  •  •  8  (  m  •  1  )  ( m  *  1 )  ■ 

Thus  if 


y(s>T>[ V(S>(1) , V<s>(2) , .  .  .,V<s>(ra)] 

and 

IT= (1(1), 1(2),.  .  . , I ( m )  | . 

The  matrix  2  can  be  decomposed  into  two  parts  as 
Z=Zq«-Zl>  where 

Zq  is  the  generalized  impedance  matrix 

and 


Z^  is  the  load  matrix. 

Assuming  that  all  dipoles  are  loaded  with  the  same 
lead  z^,  the  matrix  is  given  by 


ZL=diag(z1  7.1  .  .  .  z ! }  . 


Xi-ivi  y 2  •  •  •  ym:T- 

and 

12^1X2  X3  ■  •  '  >-(m*l)lT' 

we  can  write 


ii  - 

G11  i 

1  *  Gi2  i:2  *  ii; 

(14) 

and 

12  - 

G?  1  X 

1  ♦  Go-?  X-h  *  f«  ^  , 

(15) 

where 

Gll- 

g12  • 

G2  1  •  G22  •  tii  •  ii2' 

Xj  and  x2 

given 

by 

The  ij-th  element  of  Z  is  z  1 J = z ^ * z } 6 j j .  The 
voltages  induced  on  a  load  z^  are  given  by 

V< c>  =  ZL  I  and  I  =  ZL-J  1 > . 

however , 

v(inc)  -  zi  »  zczL-’v(0  »  y(')i 


which  implies  that 

V(t)  =[I*Z0Z l*1)'1  V(inc).  (10) 

Let  H  be  the  matrix 

H=  [ I*Z0ZL-1 ] .  (11) 

Thus,  when  incident  signals  are  impinging  on  the 
array  and  in  the  presence  of  additive  noise,  the 
output  of  the  linear  array  will  be 

V<0  -H~  1  v(  inc)  .  Ni 

For  simplicity,  let  and  K.V^).  Ue  now 

have  a  relationship  between  the  incident  signals 
and  the  received  signals  at  the  outputs  of  the  ar¬ 
ray,  which  is 

T  =  H"1  X  -  M.  (12) 

APPLICATION  TO  ESPRIT 

Consider  a  linear  array  of  (m*l)  sensors  and 
assume  there  are  d  (d<n)  narrowband  sources  lo¬ 
cated  at  angles  Oj, ;  k«l,.  .  .,d. 

First  Application 

In  the  first  application  we  consider  two 
sub-arrays  consisting  of  the  first  m  sensors  and 
the  last  m  sensors.  The  observed  signal  vector  at 
the  out;-"-  of  the  array  can  be  written  as 


Xl=(xl  x2  •  •  •  xmlT- 
-2* 1 x2  x3  •  •  •  x( m» 1 ) I  1 ■ 


On  =  (£lli  gll2  •  •  •  Slim  I. 

Isu  «2i  ••  •  Smi)1;  i-1 . m, 

Gj2T  =  (0  0  .  .  .  0  ), 

G21T  ■  (Sill  0  .  .  .  o  I, 


g22T  =  1S112  &113  '  •  •  £22(m*l)  I. 

fit! i-  Ig2i  83 i  •  •  •  8(m+ 1 ) i  1T '  u2 . 

Nj  =  | n j ,  n2,  •  .  ..  nm)T, 

i?2  -  1  n2  '  n3 . n(  m*  1 )  1T  ■ 


Consider  the  vector  Z  defined  as 
Z  =  [ Vj  Y2  )t. 

Z  can  be  written  as 


’  Gn  C12  | 

pi  j 

'  Nl  ' 

•  g2 1  g22  J 

[  x;  j 

-  iiz  ■ 

Assuming  that  the  signals  and  noise  are  statisti¬ 
cally  independent  and  that  the  noise  components 
ace  uncorrelated  from  sensor  to  sensor  with  vari¬ 
ance  a-,  Then  Czz  =  E|Z  Z1' |  is  given  by 


C 


F11 

Ui 


c12 

G22 


T(XiX,"| 

£li<li<2Hl] 

G11  g12 

E[X2XiH) 

EIX;X2»]j 

g21  g22- 

(17) 


2 


Let  (G|  and  (I)  be  the  matrices 
G11  g12 
g21  g22 


1C)  = 


and  1 1 )  = 


Im  H„ 


whera  1^  is  the  identity  matrix  and  Ilm  and  I2m 
are 


U„ 


I2„ 


0 

0 

0 

.  0 

0  • 

1 

0 

0 

.  0 

0 

0 

1 

0 

.  0 

0 

0 

0 

0 

.  0 

c 

0 

0 

0 

.  1 

0 

■  0 

1 

0 

0  . 

0  • 

0 

0 

1 

0  . 

0 

0 

0 

0 

1  . 

0 

0 

0 

0 

0  . 

1 

.  0 

0 

0 

0  . 

0  . 

Then 

1C)-1  (Czz  -a2[l ] )  ([G)-1)9 
E 1 X i  XXH ]  E l X !  X2«, 

E(X2  Xjh)  eix2  x2h] 


(18) 


Having  recovered  the  matrix  on  the  right  side  of 
equation  (18),  the  matrices  M=E!xlxlHl  and 


N=E(XiX2H|  can  be  identified.  It  can  be  shown  that 
M  and  N  have  the  decomposi t ions 

H  .  ASAh  and  N  =  AS*flAH  (19) 

where  A,  S  and  *  are  the  following  matrices 
S=E[ S  SHJ, 

ST=  1  s  i ,  .  .  •  ,  )  impinging  signal  vector, 

A  »  |ai  aj  .  .  .  aj  | 

2i  .  (1  el*i  .  . 

4>  =  diag  [  eJ  *1 ,  .  .  .,  eJ^d  ), 

♦j,  =  (wA/c)sin(e^)  ,  k = 1 ,  2,  .  .  .,  d. 
Therefore,  the  effects  of  mutual  coupling  have 


the  matrix  pencil  (M-XN)  are  given  bv 
>.  „  ej  ( wA/c )s in(  9j ) .  i  =  1(2 . d. 

Second  Application 


(20) 


In  the  second  application,  two  neighboring 
sensors  are  considered  as  a  doublet.  Assume  then 
we  have  a  linear  array  of  2m  sensors  so  as  to  con¬ 
stitute  m  doublets  and  let  there  be  d  (d<m) 
sources.  Again,  the  received  signal  at  the  output 
of  the  array  is  modeled  as 

X  «  G  X  ♦  N  (21) 
where  G  is  given  by 


'  811 

812  ’  • 

•  •  8l(2m) 

821 

822  ■  • 

•  •  82 ( 2m ) 

■  8 ( 2m )  1 

8( 2m)2  ■  • 

•  ■  8(2m)(2m)  ■ 

Let  vj  and  wj  be  the  signals  received  at  the  i-tli 
doublet .Then 


'i  -  V ( 2 i - 1 )  and  ui  ”  y<2i )  • 


(22) 


Collecting  all  the  v • ' s  in  a  vector  V  and  all  t he 
Vi's  in  a  vector  V  ,  ve  have 


V  *  Gjj  Xi  ♦  0^2  X?  ♦  Nj 

V  -  G21  xx  +  ^22  ^5  4  Nn 


and 


(??■) 

(24) 


where  G  ^ ,  Gj2*  ^21'  ^22»  X] ,  X2.  Ni  an^  *2' 
given  by 

Xl=lxl  x3  .  .  .  x(2m-l)l  ' 


X2,|x2  x4  .  .  .  x(2m)l 

G11T  =  igiil  Sll2  '  ’ 

SAAi =  lg(2i-l)l  8(21-1)3  • 
i  =  1 ,  3 . ( 2m- 1 ) , 

g12T  =  ISlil  £1^2  •  • 

£lli=  I 8f 2 i -1 )2  8(21-1)4  * 
i.l,  3 . ( 2m- 1 ) , 

T 


£  -,1  n  3  ’ 

■  8( 2 i - 1 ) ( 2m- 1 )  I  • 

£llm  I  • 

■  8(2i-l)(2n) 1 : 


•  £llm  I- 

8( 2 i ) ( 2m- 1 ) I  : 


g2  1  *  *  Ifilll  £112  • 

Slit*  1 8( 2 i )  1  8( 2i ) 3  •  • 
i  =  2,'*, .  .  . ,  (2m) , 

g22T  -  IgZh  £222  •  •  •  g22„  1. 

glli  =  l8(2i )2  8(2i)4  •  •  •  8(2i)(2m) 1 
i  =2 ,  4 . (2m), 


Nj  =  (nj,  nj, 
N2  =  l n2 ,  n4, 


n(2m-l)l 

n2mlT- 


Consider  the  vector  Z  defined  as 
Z  =  | V  U)T. 

Z  can  be  written  as 


of 

z  = 

'  Gj j  G12 

’  X1  ' 

+ 

'  Si  ' 

•  g21  g22  - 

•  x2  - 

■  N2  ■ 

(25) 


Assuming  that  the  signals  and  noise  are  statisti¬ 
cally  independent  and  that  the  noise  components 
are  uncorrelated  from  sensor  to  sensor  with 
covariance  matrix  o2I2m  where  Ijn  is  the  (2mx2m) 
identity  matrix.  Then  Czz  =  E|Z  Z™|  is  given  by 


pa 

G1 2 1 

[c2i 

g22J 

lEI^Xi”)  E|X;X2 


-  *  -  *  I 

<2X2M1- 


'll  u  1 2 

nl  Gm 


ill 


*  01 2m 

Let  [Gj  be  the  matrix 
G11  G12 


(26) 


I G ) 


Gn 


21  '••22 


3 


Then 


IGI'1  (C„  -o2I2m)  ( [G)-1  )H 
E|Xj  XjHl  ElXj  X2HI 
E|X2  xx"|  E(X2  x2h] 


(27) 


Having  recovered  the  matrix  on  the  right  side  of 
equation  (27),  the  matrices  M  =  E( X^Xj  H )  and 
N=E|X^X2^|  can  be  Identified.  Recall  that  M  and  N 
have  the  decompositions 

M  -  ASAh  and  N  =  AStHAH  ,  (28) 

where  A,  S  and  *  are  given  by 
S-=|S  §Hj 

S‘=(sj,  .  .  .  .Sj)  impinging  signal  vector, 

at  =  M 

»  =  diag  (  ej*l,  |,.  .  •,  eJ  *d  |, 

+),  =  ( w A/  c  )  s  i  n  (  0^ )  ,  k  =  1 ,  2,  .  .  .,  d. 
Therefore,  the  effects  of  mutual  coupling  have 
been  eliminated  and  the  rank  reducing  values  of 
the  matrix  pencil  (M-AN)  are  given  by 
x.  _  ej(uA/c)sin(0i).  U1(2 . d.  (29) 


COMPUTER  SIMULATION 

The  scenario  used  for  this  simulation  con¬ 
sisted  of  two  incoherent  sources  (d=2)  which  are 
incident;  on  a  linear  array  consisting  of  eight 
half  wavelength  dipoles  (m=8).  The  sources  are  as¬ 
sumed  to  be  located  at  0^=18°  and  92*22°  ,  The 
noise  was  simulated  to  be  white  Gaussian  with 
aero-mean  and  unit  variance.  The  sensors  were 
positioned  at  half  wavelength  apart  such  that 
uD/c  -  it  .  100  snapshots  were  taken  each  time  and 
the  experiments  were  repeated  50  times.  The 
results  of  the  simulation  are  shown  below, 
first  Application 


(Without  compensation  for  the  mutuals) 


SNR 

mean 

el 

mean 

02 

variance 

0i 

variance 

02 

30  dB 

13.1740 

70.4337 

1.29940 

75.60234 

25  dB 

11.9686 

36.3590 

11.87542 

31.17404 

20  dB 

14.4098 

32.0397 

6.606652 

12.23856 

15  dB 

15.4838 

31.0915 

5.512074 

11.46014 

10  dB 

15.9031 

30-8755 

7.318329 

13.04751 
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(With  compensation  for  the  mutuals) 


SNR 

mean 

0i 

mean 

$2 

vari ance 

01 

var lance 

30  dB 

18.0652 

22.2367 

0.330710 

0.  285269 

25  dB 

18.0561 

22.3543 

0.669874 

0. 748764 

20  dB 

18.2311 

22.5749 

1.683848 

1.8798H9 

15  dB 

18.3796 

22.8772 

3.496217 

4.794093 

lOdB 

18.3964 

23.5170 

8.045362 

9.578419 

Second  Application 


(Without  compensation  for  the  mutuals) 


mean 

mean 

variance 

variance 

SNR 

01 

e2 

0i 

02 

30  dB 

12.6221 

29.6850 

1.967561 

0.634513 

25  dB 

12.7642 

29. 7487 

9.552104 

3.538558 

20  dB 

14.2713 

31.4526 

35.51011 

51.34790 

15  dB 

18.4257 

40.4587 

72.28546 

242.0847 

10  dB 

19.3866 

42.4492 

66.07262 

279.2892 

(With  compensation  for  the  mutuals) 

mean 

mean 

variance 

variance 

SNR 

®i 

e2 

01 

02 

30  dB 

18.0593 

22.1633 

0.071976 

0.197370 

25  dB 

18.0902 

22.2343 

0.148269 

0.514200 

20  dB 

18.1557  J  22.4040 

0.314386 

1.135701 

15  dB 

18.2589 

22.3398 

0.737347 

2.406417 

10  dB 

17.9905 

22.8160 

2.897151 

9.306142 

Note  that  extremely  poor  estimates  are  obtained 
without  compensation  for  the  mutuals  in  both 
cases.  Compensation  results  in  significant  im¬ 
provement  . 
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